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O 
QJ ' A system consisting of two parallel coupled channels where particles in one of them 

H ■ follow the rules of totally asymmetric exclusion processes (TASEP) and in another one 

move as in symmetric simple exclusion processes (SSEP) is investigated theoretically 

Particles interact with each other via hard-core exclusion potential, and in the 

c/3 . asymmetric channel they can only hop in one direction, while on the symmetric lattice 

particles jump in both directions with equal probabilities. Inter-channel transitions 

are also allowed at every site of both lattices. Stationary state properties of the 

system are solved exactly in the limit of strong couplings between the channels. It 

is shown that strong symmetric couplings between totally asymmetric and symmetric 

C^ . channels lead to an effective partially asymmetric simple exclusion process (PASEP) 

^j \ and properties of both channels become almost identical. However, strong asymmetric 

couplings between symmetric and asymmetric channels yield an effective TASEP with 

nonzero particle flux in the asymmetric channel and zero flux on the symmetric lattice. 

J> ' For intermediate strength of couplings between the lattices a vertical cluster mean-field 

£L^ ■ method is developed. This approximate approach treats exactly particle dynamics 

qq , during the vertical transitions between the channels and it neglects the correlations 

f^*) ' along the channels. Our calculations show that in all cases there are three stationary 

r~^. . phases defined by particle dynamics at entrances, at exits or in the bulk of the 

^^ | system, while phase boundaries depend on the strength and symmetry of couplings 

between the channels. Extensive Monte Carlo computer simulations strongly support 

our theoretical predictions. Theoretical calculations and computer simulations predict 

that inter-channel couplings have a strong effect on stationary properties. It is also 

S^ ■ argued that our results might be relevant for understanding multi-particle dynamics 

^H ' of motor proteins. 
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1. Introduction 

The majority of systems in nature operate far from equilibrium, but there is no developed 
theoretical framework for comprehensive analysis of non-equilibrium processes. In this 
situation, a critical role for understanding different complex phenomena in Chemistry, 
Physics and Biology is played by a class of low- dimensional non-equilibrium multi- 
particle models known as asymmetric simple exclusion processes (ASEP) (TJ [2j [3j [U [5] . 
ASEP are stochastic models where particle interact via an exclusion potential and 
move along discrete lattices. Mechanisms of many non-equilibrium processes, such as 
biological transport, kinetics of protein synthesis and biopolymerization, car traffic, 
and hopping of quantum dots, have become better understood due to the successful 
description via asymmetric exclusion models ^l El El El EOl HH HSl H31 H31 HSl HSl [TZl HHI - 

Most theoretical studies of exclusion processes concentrate on single-lane systems 
where important exact solutions have been obtained in several cases [H H]. Recently, 
a lot of attention has also been devoted to parallel multi-chain exclusion processes 
[ISll2Dll2Ill22ll23ll23ll25ll2Sll2Zll2Sll2Sll3ni. The study of these models has been greatly 
stimulated by experimental advances in analysis of motor proteins dynamics along 
cytoskeleton filaments, in transport of mesoscopic quantum systems and in vehicular 
traffic processes [TT1 |2?1 131] - In parallel multi-chain exclusion processes particles can 
jump along the horizontal chains, but they can also switch stochastically between 
different lanes. Theoretical analysis of different multi-chain ASEP suggests that coupling 
between the channels strongly influences stationary-state phases and particle properties. 
It can produce a complex dynamic behavior, leading to many unusual phenomena, such 
localized domain walls and symmetry breaking [211 [23, [281 EH] . 

Investigations of multi-chain exclusion processes mostly involves coupling of 
asymmetric channels where the direction of particles motion is biased at each site on all 
channels. The aim of the present paper is to analyze parallel coupling of symmetric and 
asymmetric exclusion processes. The problem is motivated by the cellular transport of 
motor proteins along rigid protein filaments such as actin filaments or microtubules [31j . 
Motor protein molecules can move mostly in one direction when they are tightly bound to 
protein lattices. Occasionally, motor proteins might dissociate from the filament to the 
surrounding solution where they perform unbiased diffusional motion. Freely diffusing 
molecules can also bind to protein filaments. There is only one previous theoretical work 
that investigates coupling of TASEP and SSEP, although periodic boundary conditions 
and symmetric transition rates between the lanes are assumed [32] . Using several mean- 
field approaches and extensive computer simulations it was found that there is unequal 
redistribution of particles between different channels depending on the densities [32]. 
In a related study, Lipowsky and coworkers investigated transport of molecular motors 
in open tube that contains a single filament [HI [13]. Bound to the filament particles 
undergo asymmetric exclusion process, while the unbound molecular motors diffuse 
freely in the tube around the linear chain. The analysis of the molecular motor transport 
via tube-like compartments, performed with the help of mean-field methods and Monte 
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Carlo simulations, revealed that there are three stationary phases with phase transitions 
specified by the precise choice of boundary conditions [IT] . In our work we investigate 
a two-channel system consisting asymmetric and symmetric exclusion lanes with open 
boundary conditions, and with symmetric and asymmetric transition rates between the 
channels. 

This paper is organized as follows. In section 2 a detailed description of the 
model is given, and exact solutions for strong couplings and approximate solutions 
for intermediate couplings are presented. In section 3, we discuss Monte Carlo 
computer simulations and compare them with theoretical predictions. The final section 
4 summarizes and concludes. 

2. Theoretical analysis 

2.1. Model 

Our model consists of two parallel one- dimensional lattices as illustrated in figure 1. 
Both lanes have L sites, and we are interested in obtaining thermodynamic limit results 
when L>1. Particles move along the channels by hopping between the lattice sites 
with the exclusion potential, i.e., each site can be occupied by no more than one particle. 
In our model we apply a random sequential update when at each time step dynamics at 
randomly chosen site is followed. Particles can enter the system with the rate < a < 1 
if one or both first sites on channel 1 and 2 are not occupied. Similarly, particles exit 
the system with the rate < j3 < 1 if any of last sites are occupied. In the bulk of 
the system dynamic rules depend on the lattice: see figure 1. The particle at site i can 
switch to the same site on the lattice 2 with the rate W\ if this site is empty, or with 
the rate 1 — W\ it hops to unoccupied site i + 1 on the lattice 1. However, if the site i 
in the second channel is occupied, the particle jumps in the horizontal direction to the 
right with the rate 1 if the forward site is available. The particle at the site % on the 
lattice 2 can move vertically with the rate w 2 if the upper site i is free, or it can jump 
horizontally in either direction with the rate (1 — W2}/2 if sites % + 1 or % — 1 are available. 
However, if the upper site i is already occupied, the particle can move with the rate 
1/2 in the forward or backward directions, assuming that any of these moves are not 
blocked by already present particles at sites % — 1 or % -+- 1. The total transition rate out 
of every site i in any channel is equal to one. When the transition rates between the 
channels are equal (wi = w 2 ) the coupling is symmetric, while for w\ ^ w 2 the coupling 
between the lattices is asymmetric. 

When the inter-channel transition rates are equal to zero (u>i = w 2 = 0), the system 
decouples into two independent single-lane exclusion processes: the upper channel 
becomes a totally asymmetric process and the lower channel is a symmetric process. 
Exact solutions for single-lane TASEP and SSEP are known [U [5], and they provide 
a full description of all dynamic properties at large times. For TASEP there are three 
stationary-state phases. If the entrance to the lattice is a rate-limiting step, which 
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Figure 1. Schematic picture of two-channel system that couples symmetric and 
asymmetric exclusion processes. On the lattice 1 (upper) particles can move only 
to the right, while on the lattice 2 (lower) there is no preference in the direction of 
motion. The inter-channel transitions rates are w\ and W2- Allowed transitions are 
shown by arrows. Entrance rates at both lattices are equal to a and exit rates are 
equal to (3. 



happens for a < [3 and a < 1/2, the system is found in a low-density (LD) phase with 
the particle current and bulk density given by 

Jld = a>(l - a) , Pbuik,LD = a. (1) 

However, when exit of particles controls the overall dynamics (j3 < a and /3 < 1/2), the 
system is in a high-density (HD) phase where the stationary current and bulk density 
are 

Jhd = P(1-P), Pbuik,HD = 1-/5. (2) 

In the third phase, called a maximal-current (MC), the dynamics is governed by bulk 
processes (a > 1/2 and (3 > 1/2), and stationary properties are the following, 



J 



MC 



Tj Pbulk,MC — -• 



(3) 



The stationary properties of SSEP are much simpler with only one non-equilibrium 
phase at all conditions [5]. For the lattice with L sites the density profile is linear, 

L-i + l//3 



9 L + l/a + l/p-1' 
The average current in the steady-state of SSEP is given by 

1 



J 



(4) 



(5) 



L + l/a + l/P-1 

It can be easily seen that in the limit of L — > oo the current in the SSEP is approaching 
zero. 

There are very few exact results for multi-channel exclusion processes obtained by 
mapping them into effective single- lane exclusion models [20]. Typical approaches to 
analyze coupled multi-channel systems involve various mean-field treatments supported 
by computer simulations [191 [2Ql EH [221 HSl EH EH [26l [271 I2HI [291 E0] . Probably, one 
of the most successful approximate approaches is a vertical-cluster mean-field method 
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that describes the dynamics of inter-channel transitions explicitly and 
neglects correlations for horizontal transitions. We will also utilize this method for 
analyzing parallel coupling between TASEP and SSEP. Each vertical cluster can be 
described by introducing functions Py (i,j = 0, 1) that define the probability of different 
states. Poo corresponds to the state when both sites of the vertical cluster are empty, 
P\\ describes the state with both sites occupied, and P w and P i specify partially 
filled vertical clusters with the occupied site on the channel 1 or 2, respectively. These 
probability functions are related via a normalization condition, 

Poo + P10 + P01 + Pn = l. (6) 

2.2. Strong couplings: Exact results 

In order to understand mechanisms of coupling between TASEP and SSEP, it is 
instructive to consider strong coupling regimes. First, let us analyze the case of 
symmetric coupling with W\ = W2 = 0. In this case in the stationary-state limit it 
is not possible to observe the vertical configuration {00} with two empty sites. This is 
because this configuration can only be obtained by moving the particle horizontally if 
the previous state of the same cluster was {10} or {01}. However, the rates for these 
transitions are zero (1 — W\ = 1 — W2 — 0), and we conclude that for any bulk vertical 
cluster Poo = 0. It can be also argued that symmetry of the coupling requires to have 
P01 = P01 for any bulk vertical cluster. Then in the system there are only two types 
of vertical clusters: half- filled and fully filled. One can view {11} vertical clusters as 
effective "particles" and half-filled vertical clusters as effective "holes." These "particles" 
can advance forward with the rate p = 3/4, or they can move backward with the rate 
q = 1/4. They enter the system with the effective rate a e ff = a and leave it with the 
effective rate (3 e ff = 2/3. The factor 2 comes from the fact that in the last vertical cluster 
both particles can exit independently. Thus the two-channel system is mapped into a 
new effective single-channel partially asymmetric simple exclusion process (PASEP), for 
which exact solutions are known |33j . 

Similarly to totally asymmetric exclusion processes there are three stationary phases 
in PASEP [33]. The LD phase exists for a < (3 and a < (p — q)/2, and the stationary 
properties of this phase can be written as 

a(p-q-a) a 

■JLD — , Pbulk,LD — • [I ) 

p — q p — q 

The HD phase can be found for (3 < a and (3 < (p — q)/2 with the particle current and 
bulk density given by 

T P(p-q-P) , P (R , 

JHD — 1 Pbulk,HD — J- • \p) 

p — q p — q 

In the MC phase (a > (p — q)/2 and (3 > (p — q)/2) the stationary properties are the 
following, 

t P-Q l fa\ 

JMC — — T ) Pbulk,MC — „• {?) 
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Applying these results to our system, it is trivial to show that the LD phase is specified 
by a < 2/3 and a < 1/4 with the particle current and bulk density given by 

Jld = a(l-2a), p bu ik,LD = 1/2 + a. (10) 

The HD phase exists for a > 2/3 and j3 < 1/8 with the following stationary properties, 

J HD = 2/3(1 -4/3), Pbulk! H D = 1 - 2(3. (11) 

For a > 1/4 and (3 > 1/8 we have the MC phase with 

JhlC = 1/8, Pbulk,MC — n- (12) 

The resulting density profiles for strong coupling limit are shown in figures 2a, 3a and 
4a, while the phase diagram is outlined in figure 5a. 

Exact solutions via mapping to a single-channel exclusion process can also be found 
for strong asymmetric couplings. Let us show this for the case of Wi — 1 and W2 = 0. 
For the arbitrary bulk vertical cluster at site i we have Piq = 0, because there should 
not be the overall vertical current in the system at large times. It can be argued that the 
whole lattice 2 is fully occupied at large times after the system reaches a steady state. If 
any vacancy appears then it will be quickly filled by vertical transition form the upper 
channel. This observation suggests that bulk vertical clusters can only be found in two 
states, {11} and {01} with P 00 = P w = 0. Note, however, that other cluster states 
might exist near the boundaries, but they should not affect the overall dynamics in the 
system. We associate the fully filled vertical clusters {11} with "particles", while the 
vertical clusters {01} can be viewed as "holes." There is the particle flux in the upper 
channel, and there is no current on the lattice 2 and the bulk density in the channel 
2 is always equal to one. Thus, we mapped the system that couples asymmetric and 
symmetric exclusion process into new effective TASEP in the limit of strong asymmetric 
coupling. 

The effective entrance rate for the "particles" is equal to a e ff = a. However, the 
exit process should be considered more carefully. The overall particle flux to leave the 
channels can be written as 

Jevit = PeffPll = P(2P n + Pol) = 0(1 + Pll), (13) 

and it should be equal to the bulk current, 

J b uik = Pii(l-Pn)- (14) 

It leads to the following relation for density of fully filled vertical clusters at the end of 
the system, 

Pn = Lz^bfEMU , (15) 

which also yields the effective exit rate, 

l+/?-y//3 2 -6/3 + l 

Peff = \ • (16) 
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Now, using known results for TASEP we can predict that in the two-channel exclusion 
system with strong asymmetric coupling there are three stationary phases. Entrance 
dominated low-density phase exists for a < 1/2 and (3 > " 2 ~^J . Here the stationary 
properties are 

JlD = «(! - «)» Pb\!lk,LD = «' ( 17 ) 

Exit dominated high-density phase is found for f3 < 1/6 and (3 < a , 2 Za) w ^ n the 
following particle current and bulk density on the lattice 1: 

,(D _ 13(3-/3 + y//? 2 -6/3 + 1) (1) _ l-/? + V/?2-6/3+l 

J i?D ~~ 9 ' Pbulk,HD ~ 2 "^ 

Finally, in the maximal-current phase (for a > 1/2 and /3 > 1/6) we have 

J MC = V^, Pbulk,MC = 2' ' ' 

Density profiles for strong asymmetric coupling with wi = 1 and w 2 = are shown in 
figures 2c, 3c and 4c, and the phase diagram is presented in figure 5c. 

The other case of the asymmetric coupling, when w\ — and vj-i — 1 , can be easily 
analyzed if we recall the particle-hole symmetry of the system. The flux of particles 
moving from the left to right can be viewed as a flux of holes moving in opposite 
direction. Then stationary properties in this case can be obtained from Eqs. ( TT7|) . ( TTSj) 
and (fl9|) derived above if we perform the symmetry operations a <-> (3 and «-»■ 1. In 
this case, the nonzero particle current will be found only in the upper (asymmetric) 
channel, but there will be no particles and no flux in the bulk of the symmetric lattice. 

2.3. Intermediate couplings: Approximate theory 

When couplings between asymmetric and symmetric exclusion processes are not strong 
(w\ < 1 and/or w; 2 < 1) it is not possible to solve the system exactly via mapping 
procedure. Then an approximate theory should be developed. We will use the vertical- 
cluster mean-field approach [20], EU [29] that was successful in description of other two- 
channel exclusion processes. 

The overall properties of the system can be found by monitoring changes in four 
vertical clusters at each site. Assuming that the behavior is uniform along the lattices, 
the dynamics of bulk vertical clusters is governed by three independent master equations: 

dP n 



dt 

dPi 

dt 
dP i 

dt 



(2- Wl - w 2 )P l0 P i - 2P n P 00l (20) 

w 2 P 0l - Wl P w + 2P u Poo ~ (2 - w x - w 2 )P w P 01 , (21) 

w x P w - w 2 P m + 2P 11 P 00 -(2- Wl - w 2 )P 0l P 10 . (22) 



At large times the system reaches steady state, implying that —#■ = for i, j = 0,1. 
Then from equations (J2"T|) and (|2"2"|) we can immediately conclude that 

w 2 P 01 = WtPio. (23) 
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This expression can be understood as an equilibrium for vertical transitions between 
the channels in the bulk. Then substituting this relation along with the normalization 
condition (jHJ) into equation (120]) produces 

(2- Wl - w 2 )—P? + 2(1 + — )PnPio - 2Pn(l - Pn) = 0. (24) 

w 2 w 2 

We expect that, similarly to the cases of strong couplings, there are three stationary 

phases. To obtain stationary properties explicitly we need expressions for entrance, exit 

and bulk currents through the system, 

J entrance = «(2PoO + pLO + Pol); (25) 

■hulk = [Pi + (1 - Wi)Pio] (Poo + Poi), (26) 

J exit = (3 [2P n + (1 - Wl )P 10 + (1 - w 2 )Pqi] ■ (27) 

Substituting into these relations the values for P i an d Poo from the equilibrium for 
switching between the channels (1231) and from the normalization ([6]) we obtain 



" entrance Cx 



2(1-P n )-(1 + ^)P 10 
w 2 



(28) 



Jbuik = [Pn + (1 - wi)Pio] (1 - Pi - Po), (29) 



Jexit P 



2P n + (1 - 2 Wl + — )Po 

W2 



(30) 



The general strategy for solving the system is the following. From equation ( 1241) we 
express Pio in terms of Pn and then all stationary quantities will depend only on one 
variable. The conditions for existence and dynamic properties of entrance-dominated LD 
phase and exit-dominated HD phase can be found from the condition of the stationarity 
of the current. The MC phase can be determined by solving d ^ lk = 0. The bulk 
densities in each channel can be calculated from 

Pbuik = ^ii + ^10; Pbuik = p n + Poi- (31) 

We now proceed to analyze symmetric couplings with w\ — w 2 — w < 1. The 
equilibrium for vertical transitions gives us Poi = Pio and equation (|2"4"|) simplifies into 

(1 - w)Pl + 2P n P 10 - P n (l - P n ) = 0, (32) 

which yields the following solution, 



P 10 = V P n[^-w + wP 11 )]-P 1 ^ 
1 — w 
The properties of LD phase can be computed from the condition J en trance — Jbuik which 
leads to 



^/(l-w) 2 + lQwa 2 -{l-w) 

Pi o • l d4 J 

2w 
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The expression for the current can be written in terms of Pu as 

1 - w + wP u - y/P n (l-w + wP u 



•/ 



LD 



2a 



(35) 
1 — w 

while the bulk densities in both channels are equal to each other, and they are given by 

y/Pu[l-W + wPu) 



(1) _ (2) 
Pbulk ~ Pbulk 



wP 



11 



1 — w 

Substituting equation (134"|) into equations (135]) and (136]) we obtain 

a 



J, 



LD 



W 



and 



(1) _ (2) 
Pbulk ~ Pbulk 



y/(l - w) 2 + 16wa 2 + (l-w)-4a 
(1 - w) + 4a - y/(l -w) 2 + 16wa 2 



(36) 



(37) 



(38 



2(1 -w) 

Similarly, the properties of HD phase follow from the relation J exit = J bulk-, from 
which we obtain 

2xu(l - 2/3) - 1 + y/1 - 8wp + 16wp 2 



Pu 



1w 



The stationary current can be written as 



Jhd = 2(3y/P 11 (l-w + wP 



n , 



(39) 



(40) 



and the bulk densities are are the same as in equation ( 136|) . Using equation ( 1391) the 
explicit expressions for the current and densities are the following: 



Jhd = /? (l -4/?+ y/1 - 8w(3 + lQw(3 2 ^) , 
L-2/3. 



(1) _ (2) 
Pbulk ~ Pbulk 



(41) 
(42) 



The surprising result is that the bulk densities in HD phase, in contrast to the LD phase, 
are independent of the coupling strength. 

In the MC phase the current can be derived from equation (|29|) . 



Jmc = a/Ph(1 -w + wPn) 



1 — w + wP u 



'Puil-w + wPu 



w 



.(43) 



From the condition of maximum of the current, *? tfc = 0, we can obtain the explicit 
form of Pu for every value of the inter-channel transition rate w, 



Pi 



3w-2 64 - 144m; + lUw 2 



Qw 



+ 



+ 



12k; 



(44) 



where 



192wVY 
Y = 8 - 27w 2 + 27w 3 + 3VsVl6w - 68w 2 + ll5w 3 - 90w 4 + 27w 5 . (45) 
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Boundaries between stationary phases can be determined from the particle currents for 
each regime at transition lines. For example, HD and LD phase are separated by a curve 
given by 

a ( 1 — w — 4a + a/(1 — w ) 2 + lQwa 2 ) . . 

— ^ '- = f3 (1-4(3+ Jl-8wP+16w0 2 ) .(46) 

1 — w V / 

Note, that when w = 1 we recover the results for strong symmetric coupling obtained 

in section 2.2, as expected. 

For general case of intermediate couplings between symmetric and asymmetric 

exclusion channels (w\ ^ W2) we utilize the same approach. The solution of equation 

( }2"4"|) gives us 



-P 11 (w 1 + w 2 ) + y/( Wl + w 2 ) 2 P? 1 + 2w 2 P 11 (2w 1 - Wl w 2 - w\)(l - Puj lA _, 

wi(2 -w 2 -Wi) 
Then this equation can be used to express all stationary properties in terms of only one 
variable, Pu, and explicit calculations can be done as described above for the symmetric 
coupling. 
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Figure 2. Density profiles for the entrance-dominated (LD) phase in a two-channel 
system that couples TASEP and SSEP for a = 0.05 and (3 = 0.75. a) Wi = w 2 = 1; 
b) u>i = W2 = 1/3; c) Wi = 1 and 11)2 = 0; and d) W\ = 0.4 and W2 = 0.25. Lines 
are theoretical predictions for bulk densities, while symbols correspond to Monte Carlo 
simulations: circles represent the upper (asymmetric) channel, squares are for the lower 
(symmetric) channel. 
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Figure 3. Density profiles for the exit-dominated (HD) phase in a two-channel system 
that couples TASEP and SSEP for a = 0.75 and (3 = 0.05: a) w\ = w 2 = 1; b) 
wi = u>2 = 1/3; c) Wi = 1 and u>2 = 0; and d) w\ = 0.4 and W2 = 0.25. Lines are 
theoretical predictions for bulk densities, while symbols correspond to Monte Carlo 
simulations: circles represent the upper (asymmetric) channel, squares are for the 
lower (symmetric) channel. 



3. Monte-Carlo simulations and discussions 



We presented two theoretical approaches to investigate parallel coupling of TASEP and 
SSEP. In the first approach, the mapping of two-channel systems into effective single- 
lane exclusion models with known stationary properties has been utilized for strong 
couplings. This provides an exact description of all dynamics at large times. However, 
when the strength of the couplings between the channels was not large, we utilized 
the approximate mean-field method that neglects horizontal correlations in the system. 
In order to check the validity of the approximate method and to examine theoretical 
predictions we performed extensive Monte Carlo computer simulations. 

The obtained theoretical results are valid only in thermodynamic limit, L —> oo. 
In our simulations we used L = 1000 for each channel, although in several cases we 
checked our computations also for lattices with L = 500. It was found that computed 
dynamic properties do not depend on the size of the lattices, suggesting that finite-size 
effects are negligible in our simulations. The density profiles and the particle currents 
were calculated by averaging over trajectories that had between 2 x 10 6 and 10 8 Monte 
Carlo steps. To ensure that the system reached the stationary state, first 5% of the total 
number of steps were ignored in averaging procedures. Phase transitions between phases 
were determined by observing the abrupt changes in the density profiles for transitions 
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Figure 4. Density profiles for the maximum current (MC) phase in a two-channel 
system that couples TASEP and SSEP for a = j3 = 0.75: a) w\ = w 2 = 1; b) 
W\ = Wi = 1/3; c) W\ = 1 and w-i = 0; and d) w\ = 0.4 and wi — 0.25. Lines are 
theoretical predictions for bulk densities, while symbols correspond to Monte Carlo 
simulations: circles represent the upper (asymmetric) channel, squares are for the 
lower (symmetric) channel. 



between HD and LD phases. For boundaries between HD or LD and MC phases the 
transition points were determined by observing the saturation of the particle current. 
These procedures ensure that phase border lines are determined with precision within 
0.01 units of a and j3. 

Density profiles for different symmetric and asymmetric couplings between TASEP 
and SSEP channels are presented in figures 2, 3 and 4. In all situations excellent 
agreement between Monte Carlo computer simulations and theoretical predictions is 
observed for bulk densities. Different behavior is found for symmetric and asymmetric 
couplings between the lattices. Equal vertical transitions rates make the properties of 
both channels almost the same, with slight differences near boundaries, especially for 
LD (figures 2a and 2b) and for MC phases (figures 4a and 4b). Increasing the strength 
of the symmetric coupling puts more particle in the system, and even for the LD phase 
the bulk densities are larger than 1/2: see figures 2a and 4a. Surprising results are 
found for the HD phase where bulk densities are functions of only the exit rate f3 and 
they are independent of the strength of the coupling. For asymmetric coupling densities 
in the upper and lower channels differ significantly. When the vertical transition rate 
from the TASEP to SSEP is larger it leads to larger densities in the symmetric lattice, 
and in the strong coupling limit (w\ = 1 and u>2 = 0) it even fills the second channel 
completely in all phases. For intermediate asymmetric couplings both channels behave 
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Figure 5. Phase diagrams of a two-channel system that couples TASEP and SSEP: 
a) w\ = W2 — 1; b) W\ — w 2 = 1/3; c) Wi = 1 and w 2 = 0; and d) W\ = 0.4 
and W2 = 0.25. Symbols correspond to Monte Carlo simulations, while theoretical 
predictions are represented by lines. 



qualitatively similarly. 

Phase diagrams for a two-channel systems that couples TASEP and SSEP are 
illustrated in figure 5. In all cases the system can exist in one of three stationary 
phases: the entrance-dominated low-density phase, the exit-dominated high-density 
phase and the maximal-current phase specified by bulk dynamics. Comparison between 
computer simulations and theoretical calculations suggests that our theoretical method 
quantitatively correct in description of stationary properties of this system. However, 
there are several small deviations between theoretical predictions and Monte Carlo 
results, especially for LD/MC phase transitions line for intermediate couplings (see 
figure 5d), indicating that correlations inside the lattice are important for some ranges 
of parameters. It can be seen that symmetric couplings between the channels decrease 
the phase volume for the high-density phases, while asymmetric couplings have the 
same effect on the LD phases. In the case of unequal vertical transition rates and for 
intermediate symmetric couplings the boundaries between the LD and HD phases are 
not linear, as found for strong symmetric couplings, but rather slightly curved. 

Inter-channel particle transitions influence the overall current through the system 
as shown in figure 6 for the MC phases. For any symmetric couplings the particle 
fluxes through the system will go down, as was found before in the case of two-channel 
TASEP systems [20]. However, breaking the symmetry in the vertical transition rates 
actually leads to increase in the particle current. As illustrated in figure 6, for W\ — 1 
lowering the transition rate w<i from 1 to zero increase the particle current in two 
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times. Similar behavior is observed for other stationary phases. Since the two-channel 
system that couples asymmetric and symmetric exclusion processes might be relevant 
for understanding transport of motor proteins along protein filaments [H] , we can argue 
that this observations might be important for understanding motor protein's dynamics. 
One can suggests that the flux of molecular motors can be controlled by modifying 
the association and dissociation rates to protein filaments, e.g., via changing the ionic 
strength, temperature or viscosity. 



J 0.2 




Figure 6. Particle current through a two-channel system that couples TASEP and 
SSEP as a function of the vertical transition rate wi for the fixed vertical transition 
rate w\ = 1 and for a = P = 0.75. Symbols are from Monte Carlo simulations and 
lines arc the result of theoretical calculations. 



4. Summary and conclusions 



We investigated the two-channel system that couples asymmetric and symmetric 
exclusion processes for different inter-channel transition rates in steady-state regime. 
In the limit of strong symmetric coupling, w± = u>2 = 1, the exact description of 
particle dynamics is achieved by mapping the two-channel system, consisting of totally 
asymmetric and symmetric exclusion processes, into an effective one-channel partially 
asymmetric exclusion process with known stationary properties. Exact solutions are 
also obtained in the limit of strong asymmetric inter-channel rates (w\ = 1 and u>2 = 0, 
or W\ = and iw 2 = 1)- In this case the two-channel system mapped into a single-lane 
totally asymmetric process with explicit description of all dynamic properties. 

The two-channel system with coupled TASEP and SSEP for intermediate vertical 
transition rates has been analyzed via an approximate theoretical approach. In this 
method the dynamics of vertical inter-channel transitions is fully accounted, while the 
correlations along the horizontal lattices are neglected. The vertical-cluster mean-field 
approach allowed us to calculate analytically or numerically exactly particle currents, 
bulk densities and phase diagrams. The predictions of the approximate method are 
in excellent agreement with extensive Monte Carlo computer simulations. Theoretical 
calculations and computer simulations indicate that the strength and symmetry in the 
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vertical transition rates have a strong effect on the overall particle dynamics. Equal 
inter-channel transition rates symmetrize the particle properties in both channels. 
Symmetric couplings also lower the particle fluxes and increase the bulk densities. 
Asymmetric inter-channel transitions generally lead to similar qualitative properties in 
both channels, although with different values of bulk densities and currents. Asymmetric 
couplings also increase the transport capability of the system. 

We discussed the relevance of the results of this investigation for understanding 
mechanisms of motor protein's motion along the protein filaments. It can be concluded 
from our theoretical analysis that dynamics of motor proteins can be controlled and 
modified by changing the association and dissociation rates. It will be interesting to 
test our predictions in experimental studies. 
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